LIMITING SPECTRAL DISTRIBUTION FOR WIGNER 
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^. Abstract. In this article we show the existence of limiting spectral 

C3 ' distribution of a symmetric random matrix whose entries come from a 

stationary Gaussian process with covariances satisfying a summability 
condition. We provide an explicit description of the moments of the 
£/") ' limiting measure. We also show that in some special cases the Gaussian 

assumption can be relaxed. The description of the limiting measure can 
also be made via its Stieltjes transform which is characterized as the 
_> ' solution of a functional equation. In two special cases, we get a descrip- 

J™^ , tion of the limiting measure - one as a free product convolution of two 

I" H ■ distributions, and the other one as a dilation of the Wigner semicircular 

^ C"| law. 

-(— > ' 






pq ■ 1. Introduction 

> ■ 
t-i- . In his seminal paper, Wigner (1958) showed that for a symmetric ran- 

C^S \ dom matrix with independent on and off diagonal entries satisfying some 

moment conditions, the empirical spectral distribution (henceforth ESD) 
converges to the Wigner semicircle law (defined in (7.2), henceforth WSL). 
^T ' Subsequent work has tried to obtain a better understanding of the spec- 

~f~, . trum of such matrices, which plays an important role in physics as well as 

other branches of mathematics such as operator algebras. Recently, there 
has been interest in how far the independence assumption and the moment 
conditions can be relaxed. The reader may refer to the recent review ar- 
ticle by Ben Arous and Guionnet (2011) and the references therein for an 
C^ ■ overview of currently available results. 

Relaxation of the independence assumption has been investigated by 
Chatterjee (2006), Gotze and Tikhomirov (2005), Hofmann-Credner and Stolz 
(2008), Rashidi Far et al. (2008). Adamczak (2011), Pfaffel and Schlemm 
(2012), and Hachem et al. (2005) have studied the sample covariance matrix 
imposing some dependence on the rows and columns. However, the limiting 
spectral distributions (henceforth LSD) obtained by considering symmetric 
matrices with the independence assumption weakened have stayed within 
the WSL regime for the most part. One exception is Anderson and Zeitouni 
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(2008), who considered the LSD of Wigner matrices where on and off diag- 
onal elements form a finite-range dependent random field; in particular, the 
entries are assumed to be independent beyond a finite range, and within the 
finite range the correlation structure is given by a kernel function. 

Motivation. We begin with a few examples to motivate the problem stud- 
ied in this article. In each of the following examples, a random field {Zij : 
i,j > 1} is developed. For n > 1, let A n be the n x n matrix whose (i,j)-ih 
entry is Z iA j iiV j. The question is whether the ESD of A n / y/n converges as 
n—>oo, and if so, where. 

Example 1. Let {Zij : i,j > 1} be a mean zero Gaussian process such that 
E[ZijZ i+ k t j + i] = pl fe l+rl for integers i,j,k,l such that i,j,i + k,j + I > 1, 
where \p\ < 1 is fixed. This process can be thought of as a "two dimensional 
AR(1) process", because {Z^j : j > 1} is an AR(1) process for fixed i, as is 
{Zij '■ i > 1} for fixed j. 

Example 2. Assume that {Gij : i,j > 1} are i.i.d. standard Gaussian 
random variables, and N is a fixed positive integer. Define 

N N 



■^m :— z_^ Z-/ Gi + kj+i, i,j>l- 



fc=0 1=0 



Example 3. Suppose that (G n : n € Z) is a mean zero variance one 
stationary Gaussian process. Let (G l n : n G Z) be i.i.d. copies of (G n : n G Z) 
for i = ...,-2,-1,0,1,2,.... Set Z*j := G* _i , i,j G Z. 
Example 4. Let {c^^} be real numbers such that 



5Z X] C M < °°' 



-oo l=— oo 

Cfc ; = Q t for all k, I G 



oo 



^ Cfc.jCfe/^ = for all k ^ k' . 

/ = — oo 

As in Example 2, let {Gjj '■ i,j > 1} be i.i.d. standard Gaussian random 
variables. Define 

^hi := Z_^ Ck,iGi-k,j-i, i,j G Z. 

It is shown later in Section 7 that for Examples 1 and 2, the LSD of A n /^/n 
is the free product convolution of the WSL with a distribution supported on a 
compact subset of [0, oo), and for Examples 3 (under the additional assump- 
tion that X^^=i \E(GoG n )\ < oo) and 4, the LSD is a dilation of the WSL. 
To the best of our understanding, Example 2 is the only one of the above 
examples where the result follows from the work of Anderson and Zeitouni 
(2008), because that is the only example where two entries are independent 
if their distance is above a threshold. 
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Motivated by these examples, this article considers a random matrix 
model where on and off diagonal entries form a stationary Gaussian field, 
with the covariance of the entries satisfying a summability condition. It 
is shown using the method of moments that the ESD converges to a non- 
degenerate measure. The combinatorial approach we have adopted for cal- 
culating the traces of powers of the matrices avoids the use of independence 
in any stage. Unlike Anderson and Zeitouni (2008) where the use of inde- 
pendence facilitates the negligibility of certain partitions, sharper estimates 
are needed on a class of partitions. These sharper estimates on the set of 
partitions and Wick's formula are used to derive the limiting moments. The 
assumption of Gaussianity, although important in the proof, is relaxed to 
allow for a fairly general class of input sequences using the Lindeberg type 
argument developed in Chatterjee (2005). An interpretation of the limit- 
ing moments in terms of functions of non-crossing pair partitions is used 
to derive the Stieltjes transform of the measure. The form of the Stieltjes 
transform indicates a relationship with operator-valued semicircular vari- 
ables studied in Speicher (1998) (for the application of free probability to 
random matrices, see the recent review by Speicher (2011)). 

Outline of our contribution. Let (Zij : i,j G Z) be a stationary, mean 
zero, variance one Gaussian process. Stationarity here means that for k, I G 
Z, 

(Z i+ktj+ i : i, j G Z) = (Zij : i, j G Z) . 

For i,j > 1, set 

and let 

(1-1) A n := ((X id )) nxn , n > 1 . 

Let Ai < . . . < A n be the eigenvalues of A n , which are real because A n is 
symmetric, and denote 



1 " 
(!-2) ^ : =-E^ 






The main result of this article is Theorem 2.1, stated in Section 2 along with 
an outline of the proof, which gives a set of conditions on the covariance of 
{Xij} under which the ESD [x n converges weakly in probability. In Section 
3, some combinatorial results are proven, which are used in Section 4 for 
the proof of Theorem 2.1. In Section 5, we show that by specializing on an 
infinite order moving average process with independent inputs satisfying the 
Pastur condition, Theorem 2.1 and an invariance principle can be used to 
establish the convergence of the ESD. In Section 6, an explicit description 
of the Stieltjes transform is provided using the moment formula and some 
properties of the Kreweras complement. In Section 7, two explicit examples 
are described where we get better descriptions of the limit: Theorem 7.1 
gives conditions under which the LSD is the free multiplicative convolution 
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of the WSL and another distribution. Theorem 7.2 gives conditions under 
which the LSD is the WSL. Finally, in Section 8, Theorem 7.1 is extended 
to entail the case where the correlations are not necessarily summable. That 
assumption is replaced there by the weaker assumption of absolute continuity 
of the spectral measure. 



2. The main result 

In this section, we state the main result, and give an outline of the proof. 
Let the n x n random symmetric matrix A n be as in (1.1), and set /i n to be 
ESD of A n j '\/n, as defined in (1.2). Before stating the main result, we need 
a few more notations and assumptions. Define 

(2.1) R(u, v) := E [Zo,oZ_ U)V ] ,u,v£Z. 

The assumptions are the following. 
Assumption 1: R(-,-) is symmetric, that is, 

(2.2) R(u, v) = R(v, u) for all it, v G Z . 
Assumption 2: R(-, •) is absolutely summable, that is, 

(2.3) R:= Y^ \R( u ,v)\ <oo. 

An immediate consequence of Assumption 1 and stationarity is that 

(2.4) R(u,v) = R(-v,-u),u,veZ. 

A consequence of Assumption 2 is the following. Fix a £ NC2(2m), 
the set of non-crossing pair partitions of {1, . . . , 2m}. Let (Vi, . . . , V m+ \) 
denote the Kreweras complement of a, which is the maximal partition a of 
{1, . . . , 2m} such that cU a is a non-crossing partition of {1, 1, . . . , 2m,, 2m}. 
For 1 < i < m + 1 , denote 



(2.5) 




Vi:={v\,. ..,<}. 




Define 






(2.6) S(a) := < 


(h, 


Is 

. . . ,k 2m ) e Z 2m :^2k v s = 0, s = 1, . . 
3=1 


. ,m + 1 


If a = {(ui,u m+ i),... 


, (u m ,U2m)}, then notice that 








(k 1 ,...,k 2m )eS(cr) {u,v)&a 
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E 



#a£ S(a) : there exists a permutation tt of {1, ... , m} 



such that k Uj = i^) and k Um+j 



-7t(J)' ^ 



l,...,m| 






< 



E 



m! J J |i?(zj,i m+ j; 



(2.7): 

In view of the above calculation, it makes sense to define 

(2.8) /3 2m := J2 Yl II R(k u ,h),m>l. 

<rGNC 2 (2m) keS{a) (u,v)€cr 

The main result of this article is the following. 

Theorem 2.1. Under Assumption 1 and Assumption 2, \x n converges weakly 
in probability to a distribution fi. The k-th moment of [x is zero if k is 
odd, and /3k if k is even. Furthermore, [i is uniquely determined by its 
moments, that is, if a distribution has the same moments as that of \x, then 
the distribution equals \i. 

Remark 1. As is common in the literature, the phrase "fi n converges weakly 
in probability to a distribution fi " means that 

P 

L(/i n ,/i) — > 0, 

as n —> oo, where L, the Levy distance, is defined by 

(2.9) L(y\,V2) := inf< e > : v\ ((— oo,x — e]) — e < v 2 ((— oo,x]) < 

v\ ((— oo, x + e]) + e for all x £ R > , 
for probability measures v\, v 2 on R. 

We end this section with a brief outline of the proof of the above result. 
As is standard in a proof by the method of moments, what needs to be 
shown is that for fixed m > 1, 

n 

(2.10) lirn n~( m+1 ) £ E [X h , 2 . . . X i2m _ lMm X i2mM ] = /3 2m . 

U,— ,J2m=l 

As in the proof of the classical Wigner's result, the first step is to get rid 
of the "non-pair matched" tuples i = (ii, ■ ■ -i^m) i n the above sum. Fix 
N > 1, and say that a tuple i S {1, . . . ,n} 2m is A^-pair matched if there 
exists a pair partition -/r of {1, . . . , 2m} such that for all (u, v) £ tv, 

\i u -i A i u - iv-x A i v \ V \i u -x Vi u - i v -\ V i v \ < N , 
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with the convention iq := i<i m - It needs to be shown that if Cn,u denotes 
the set of tuples in {1, . . . , n} 2m which are not TV-pair matched, then 



(2.11) lim limsupn" (m+1) 



iSCjv, 



E[X h 



i 2 . . . Aj 2r] 



,x. 



l,«2m^ i *2m,*l 



0. 



Unlike in the classical Wigner's result, this is a non-trivial step in our sit- 
uation because not only does the above sum not vanish for iV large, even 
showing that the expectation in modulus is less than some e is not enough 
because #CV in ~ n2m as n — > oo, and the sum is scaled only by n^ m+l >. 
This is precisely the step where Assumption 2 plays an important role. 

Once (2.11) is established, what remains to be shown for (2.10) is that 
for fixed N, 



(2.12) 



lim n~( m+1 ) 



iecir. 



E \_Xi t ,i 2 ■ ■ ■ X 



'2)1 



-li'2m^»2m,'l 



e e n R ^u,k v ). 

o-£iVC2(2m) fceS , (cr):maxj \kj\<N (u,v)Ecr 

By standard combinatorial arguments, the sum over C^ can be shown to 
be asymptotically equivalent to the sum over all tuples that are Catalan 
with respect to some a £ NC2(2m), that is, whenever (J, k) G a, 

\ij-l -ik\ V \ij -ifc-i| < N . 

The final step is to show that for fixed a £ NC2(2m), if D a denotes the set 
of tuples in {1, . . . , n} 2m which are Catalan with respect to <r, then 



lim J2 E i X h,i2---Xi: 






Xi 



2m,«l. 



e n R ( k u,k v )- 

keS(a):ma,Xj\kj\<N (u,v)£cr 

This follows by computing the expectation via Wick's formula, and observing 
that in that formula, the contribution of all the pair partitions excluding a 
is asymptotically negligible. This final step establishes (2.12). 

3. Some combinatorics 

In this section, we recall some elementary combinatorial notions, and 
prove a few results related to them. The results of this section are not of 
independent interest, but will be used in Section 4. 

There is an infinite totally ordered set called the "alphabet" whose ele- 
ments are called "letters" . The order with which the alphabet is endowed is 
the "alphabetical ordering" . A "word" is an ordered finite collection of not 
necessarily distinct letters. While the actual description of the alphabet is 
irrelevant, to fix ideas, we shall consider the set of natural numbers with the 
natural ordering to be the alphabet. Two words are "distinct" if one cannot 
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be obtained from the other by relabeling letters. For example, the words 
1112 and 2224 are not distinct, and the words 4612 and 2181 are distinct. 
If the lengths of two words are different, then they are necessarily distinct. 
Let P{2m) and NC2{2m) denote the set of pair partitions and non- 
crossing pair partitions of {l,...,2m} respectively. Clearly, NC2(2m) is 
a proper subset of P(2m) for m > 2. For example, {(1,3), (2, 4), (5, 6)} G 
P(6) \ NC 2 (6) and {(1,4), (2, 3), (5, 6)} G NC 2 (6). Recall that for a G 
NC2(2m), the Kreweras complement K(a) is the maximal partition a of 
{1, . . . , 2m} such that cU a is a non-crossing partition of {1, 1, ... , 2m, 2m}. 
For example, 



K ({(1,4), (2, 3), (5, 6)}) ={(1,3), (2), (4, 6), (5)}. 

For m > 1, 7r is a "pairing" of {1, . . . , 2m} if it is a permutation of that 
set satisfying 

7r(i) / i = 7t(t(j))) 1 < i < 2?n . 
Call -/r to be an "almost pairing" of {1, ... , 2m} if it is a permutation satis- 
fying 

#{l<j<2m:vr(i)=i}=2, 
and 

TrMi)) = J, 1 < J < 2m. 
There is a clear bijection between the set of pairings and P(2m), namely for 
any pairing ir, {(i, vr(j)) : 1 < j < 2m} G P(2m). Keeping this bijection in 
mind, we shall use the words pairing and pair partition interchangeably. 

Recall that a word A := a\ . . . a,2 m is "pair matched" if there exists a 
pairing it of {1, ... , 2m} such that 

aj = a„(j), 1 < j <2m. 

If there exists an almost pairing it satisfying the above, then A is "almost 
pair matched". Examples of pair matched words are 1212, 122221 etc., and 
that of almost pair matched words are 1231, 2111, 1221 etc. An example of 
a non-pair matched word is 111222. Notice that while the set of pairings and 
almost pairings are disjoint, a pair matched word is necessarily almost pair 
matched. Another useful observation is that a word is almost pair matched 
if and only if it be can made pair matched by changing at most one letter. 

Recall that a pair matched word is a "Catalan word" if successive deletions 
of double letters lead to the empty word. Examples of Catalan words are 
1221, 123321, 122122 etc., while 1212 is an example of a pair-matched word 
which is not a Catalan word. 

The conventions that we now discuss will be assumed throughout the 
article. Any tuple i G Z, k is taken to be of the form i := {i\, . . . ,ik), and 
furthermore, implicitly defines i$ := i}~. The same convention also applies 
to words, that is, for a word A = a\ . . . a^, ao '■= a>k- The next convention is 
that, for alphabets or integers a, b, c, d, we say 

(3.1) (o,6)«(c,d) 
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ifaAb = cAd and a V b = c V d. For integers "A" and "V" have the usual 
interpretation of minimum and maximum respectively, whereas for letters 
u and v, u A v and u\/ v mean the one that comes first in the alphabetical 
ordering, and the one that comes later, respectively. 

Given words A := a± . . . a m and B := b\ . . . b m of the same length, say 
that B is an "offspring" of A if the following is true. Whenever clj = a k for 
some 1 < j, k < m, it holds that 

(bj-t,^) « (6fc_i,6fc) • 

For example, 1213 is an offspring of both 1221 and 5789, and both 7565 and 
1111 are offsprings of 2211. 

The next notion we need is that of a "compound offspring word" . A word 
B = b\ . . . &2m is a compound offspring of A = a\ . . . aim if b\ . . . b m and 
b m+ i . . . bim are offsprings of a\ . . . a m and a m+ i . . . a2 m respectively, and 
furthermore, whenever cij = ak, it holds that 



( b J> b iU)) ~ i b k,b l{k) ) 



where 



j-1, jG{l,...,2m}\{l,m + l}, 
l{j) ■= { m, j = 1, 

2m, j = m + 1 . 

For a word A, we denote by ^A the number of distinct letters in A (and 
not the length of ^4). 

The following result is well known in the literature, but in different set- 
tings. One can look at, for example, equation (34) in the proof of Theorem 4 
in Bose and Sen (2008) where the same claim has been restated in a slightly 
different language. Hence we omit the proof. 

Lemma 3.1. (a) Let A be a pair matched word of length 2m for some 
m > 1. Then A has an offspring word B with 

(3.2) #B = m + l 

if and only if A is a Catalan word with 

jfA = m. 

In this case, the offspring word is unique upto relabeling of letters. 

(b) Assume that A\ and Ai are distinct Catalan words of length 2m, that 
is, one cannot be obtained from the other by relabeling letters. If B\ and B2 
are offsprings of A\ and A<i respectively such that 

jfB l = #B 2 = m+l, 

then B\ and B<i are distinct. 

(c) Furthermore, if B = b\ . . . 62m is an offspring of A = a\ . . . aim satisfying 
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(3.2), then it is necessary that whenever clj = a k for some j < k, it holds 
that 

(3.3) bj-i = bk , and 

(3.4) bj = b^. 

As mentioned at the beginning of the section, the remaining results will 
be used for later results in this section or the ones in Section 4. 

Lemma 3.2. Let B = b\ . . . b 2 k be a compound offspring word of A = 
a i ■ ■ ■ a 2k f or some k > 2. Define 

#(b k h . . . bj) - #(ai ... a,), 1 < j < k - 1 , 
#(b 2k h . . . bj) - #(oi . . . Oj), k + 1 < j < 2k - 1 . 

Then, 

(3.5) 1 > m > . . . > n fc _i > n k+1 - 1 > . . . > n 2k -i - 1 • 

Consequently, if D is an offspring word of C := C\. . .c k for some k > 2, 
then 

(3.6) #D<l + #(ci...c fc _i). 

Proof. It is easy to see that the inequality n k ^\ < 1 in (3.5) implies (3.6). 
So the former claim is the one that needs a proof. The leftmost inequality 
in (3.5) is trivial. For the subsequent inequalities, fix 1 < j < k — 2, and we 
shall show that 

(3.7) n j+i ^ n j ■ 
The proof will be separate for the two cases: 

(Case 1) Oj+i 7^ ai for all 1 < i < j , 

and 

(Case 2) aj+\ = a^ for some 1 < i < j . 
Observe that in Case 1, 

rij+i = #(b k bi . . . bj+i) - #(oi . . . Oj + i) 

< 1 + #(&fe6i • • ■ bj) - #(oi • • • o i+ i) 

= #(b k b 1 ...b j )-#(ai...a j ) 

= n j- 

In Case 2, if z > 2, then bj + \ equals bi or 6j_i, and if i = 1, then bj+i equals 
bi or b k . Hence 

#(b k b 1 ...b j+x ) = #(b k h... bj), 

and #(ai . . . a j+ i) = #(ai . . . a-,) , 

which shows that 

n j+1 = nj . 
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This establishes (3.7) for 1 < j < k — 2. Similar arguments establish the 
same claim for k + 1 < j <2k — 2, and that 

rik+i < 1 + n k -i . 
This completes the proof of (3.5), and thereby establishes the lemma. □ 

Lemma 3.3. Suppose that A is an almost pair matched word of length 2m 
with m>2, and B is an offspring of A. Then, 

#B < m + 1 . 

Proof. Since A is almost pair matched, clearly j^A < m + 1. If #A < m, 
then by (3.6), it follows that 

#B <l + #A<l + m. 

So, without loss of generality, let us assume that j^A = m + 1. We 
start with the observation that if b\ . . . b 2m is an offspring of a± . . . a2 m , then 
hi . . . h2 m h\ is an offspring of 02 . . . a2mOi- Therefore, once again without 
loss of generality, we can and do assume that 

A = W lC W 2 d, 

where W\ and Wi are possibly empty words such that W\W 2 is pair matched, 
and c, d are distinct letters which do not occur in WiW 2 - Therefore, by (3.6), 

#B < l + #(W lC W 2 ) 
= 1 + m . 

This completes the proof. □ 

Lemma 3.4. Suppose that A is an almost pair matched word of length 2m 
with m>2, and B is a compound offspring of A. Then, 

#B < m + 2 . 

Proof. Denote A = a± . . .a 2m . By the fact that the rightmost quantity in 
(3.5) is at most 1, it follows that 

(3.8) #£<2 + #(ai...a 2m _i). 

Thus, the claim of the lemma follows if #A < m. So assume without loss 
of generality that 

#A = m + 1. 

It is easy to see that since h\ . . . b 2m is a compound offspring of a\ . . . a 2m , 
so are b m+ i . . . b 2m b\ ...b m and h . . . b m b m+2 . . . b 2m b m+ i of 
a m+ i . . . a 2m ai ...a m and ai . . . a m a m+2 . . . a 2m a m+ i respectively. There- 
fore, as in the proof of Lemma 3.3, we assume without loss of generality 
that 

(3.9) a 2m j=- a,j, j = 1, . . . , 2m - 1 . 
Clearly, (3.8) and the observation that under this assumption 

#(ai ...a 2m _i) =m, 
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completes the proof. □ 

Lemma 3.5. Suppose that A is a pair matched word of length Am for some 
m > 1 . Then, A has a compound offspring word B with 

(3.10) #£ = 2m + 2 
if and only if 

(3.11) A = A X A 2 

where A\ and A 2 are Catalan words of length 2m with no common letters, 
and 

(3.12) #A 1 = #A 2 = m. 

In this case, 

B = B\B 2 , 

where B\ and B 2 are offspring words of A± and A 2 respectively, do not have 
a common letter, and satisfy 

(3.13) #J3 1 = #B 2 = m + l. 

Furthermore, the compound offspring word B satisfying (3.10) is unique up 
to relabeling, and if B and B' are compound offspring words of distinct pair 
matched words of length 4m satisfying 

#B = #B' = 2m + 2 , 

then B and B' are distinct. 

Proof. Assume for a moment that the "if and only if" claim has been shown. 
Let B be a compound offspring of A such that (3.10) holds. By definition 
of a compound offspring word, we can write 

B = B\B 2 

where B\ and B 2 are offspring words of A\ and A 2 respectively. Now notice 
that by (3.12) and Lemma 3.4, 

#Si<m + l, i = l,2. 

Thus, 

2m + 2 = #B 

< #Bi + #B 2 

< 2m + 2 . 

Therefore, B± and B 2 cannot have a common letter, because otherwise, the 
inequality in the second line becomes strict. Also, the inequality in the last 
line must be an equality, proving (3.13). The final claim follows from Lemma 
3.1 (b). 

So the "if and only if" claim is the only part that needs a proof. Once 
again, the "if" part follows trivially from Lemma 3.1 (a). Let us proceed 
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towards the "only if" part. So assume that A has a compound offspring 
word B such that (3.10) holds. Let A := a\ . . . a^ m and B = b\ . . . b^ m - Set 

A\ := a\ . . . a 2m , 

A2 := 02m+l • • .04 



in ■ 



Thus, (3.11) trivially holds. We start with showing that A\ and A 2 do not 
have a common letter. Assume for the sake of contradiction that they have 
a common letter. The arguments that justify the assumption (3.9) in the 
proof of Lemma 3.4 show that in this case, a2 m +i can be chosen to be that 
letter without loss of generality, that is, 

(3.14) o,2m+i = clj for some 1 < j < 1m . 
Therefore, 

(3.15) (6 4m , b 2m+1 ) ~ (bj,b^Q)) , 

where "«" and 7(-) are as in (3.1) and the definition of a compound offspring 
word respectively. Let rij for j = 1, . . . , 2m — 1, 2m + 1, . . . , Am — 1, be as 
in the statement of Lemma 3.2 with k = 2m. By that result, it follows that 



and 
Thus, 



1 > m > . . . > n 2m -i , 

"-2-m+l > • • • > "<4m-l • 

"-4m- 1 < "2m+l 

= #(bimbl ■ ■ ■ &2m+l) — #( a l • • • a 2m+l) 
= #(&1 • ■ ■ b2m) - #(ai ■ ■ ■ 02m+l) 

< #(b 1 ...b 2m ) -#(a 1 ...a 2m -i) 

= "2m-l 

< 1, 

the equality in the third line following by (3.15). However, notice that 

(3.16) n 4m _i = #B- #(01 . . . a 4m _i) = #5 - #A > 2 , 

the second equality following from the fact that A is pair matched, and the 
inequality following from (3.10). This clearly is a contradiction, thus showing 
that A\ and A 2 have no common letters. An immediate consequence is that 
A\ and A 2 are pair matched words. 

Next, we proceed towards showing (3.12). Lemma 3.2 implies that 

(3.17) 1 > n 2m _i > n 2m+1 - 1 > n 4m _i - 1 > 1 , 

the rightmost one following from (3.16) which is clearly valid regardless of 
the assumption (3.14). Thus, 

#A = #B-rn m -i =2m. 
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Hence each letter in A conies exactly twice, and so (3.12) holds. Another 
consequence of (3.17) is that 

nim-\ = 1 , 
a restatement of which in view of the fact that A\ is pair matched is 

#(6i • • • b 2m ) = 1 + #Ai = m + 1 . 

Since b\ . . . 62m is an offspring word of Ai, by Lemma 3.1, it follows that 
A\ is a Catalan word. A similar argument holds for A2, and completes the 
proof. □ 

Lemma 3.6. Let A := A\ . . . A2 m be an almost pair matched word of length 
2m where m > 2. Assume that A has m + 1 distinct letters a±, . . . , a m +i 
with each 0/01, . . . , a m _i occurring twice. Fix 

u := (ui,...,Um-i),v := (v±, . . . ,u m -l) G Z m_1 . 

Given a 2m-tuple (i±, . . . , iim) i n N 2m , say that it is uv-matched if the fol- 
lowing is true: 

whenever Aj = A^ = ai for some 1 < j < k < 2m and 1 < I < m — 1 , 

(ij_i A ij) - (i k -i A i fc ) = ui , 
and 

(ij_i V ij) - (i k -i V i k ) = vi , 
where io := i2m; os usual. Let U n denote the set of uv-matched tuples in 
{l,...,n} 2m forn> 1. Then, 

#U n < 4 m n m+1 for all n > 1 . 

Proof. Let IT denote the set of all functions from {1, 2, . . . , 2?n} to {0, 1}. 
Clearly, if (i±, . . . , i2m) G C/ n , then the following is true. There exists it E Id 
such that whenever A,- = A k = ai for some 1 < j < /c < 2m and 1 < / < 
m — 1, it holds that 

(3.18) *i-7i-(j) ~ ik-ir(k) = u i . 

and 

(3-19) ^■_(i_ 7r ( J -)) - ifc_(i_fl-(jfc)) = vi . 

It is easy to see that the number of (ii, . . . , i-im) i n {1> ■ ■ ■ ,n} 2m satisfy- 
ing (3.18) and (3.19) is at most the number of those satisfying the same 
equations with u\ and v\ replaced by for all I. 

Fix 7r G n. We shall now show that the number of i := {i\, . . . , i2m) G 
{l,...,n} 2m satisfying (3.18) and (3.19) with ui and v\ replaced by for 
all / is at most n m+1 . Clearly, for any such i, the word B = i\i2 ■ ■ ■ 12m is 
an offspring word of A. Since B can have at most m + 1 distinct letters by 
Lemma 3.3, it follows that the number of such i's is at most n m+1 . Thus, 

#U n < n m+1 #II = 4 m n m+1 , 
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which completes the proof. □ 

Define a binary operation i on Z 2 , that is, a function from Z 2 x Z 2 to Z 2 
as follows: 

(3.20) (i,j)*(k,l) := (iAj -kAl,kV 1-iV j), i,j,k,l G Z. 

Fix 

u := (ui,...,u m -i),v := («i, . . . ,f m -i) G Z m_1 . 

For n > 1, let V^(it, v) denote the set of all tuples i in {1, ... , n} 2m for which 
there exists an almost pairing tt of {1, . . . , 2m} and an onto function <j> from 
iy := {1 < j < 2m : ir(j) / j} to {1, . . . , m - 1} such that for all j G W, 

(3.21) <j>(j) = </>(*(j)) , 
and 

(3-22) (*j-i,ij) * (V(j)-i' *«■(?•)) = (^0>^(i)) • 

Lemma 3.7. There exists a finite constant C(m) depending only on m such 
that 

(3.23) #KK«) < C(m)n m+1 , n > 1 . 

Proof. Since there are only finitely many almost pair matched words of 
length 2m and each of them has finitely many offspring words, there are 
only finitely many almost pairings n, and given any tt the number of func- 
tions 4> satisfying (3.21) and (3.22) is also finite, the proof follows from the 
conclusion of Lemma 3.6. □ 

4. Proof of the main result 

For the proof of Theorem 2.1, we shall need a few notations. Fix JVgN. 
Call a 2m-tuple i := (i\, . . . , i^m) G N 2m "iV-Catalan corresponding to cr" if 
there exists a G NC2(2m) such that whenever (j, k) G o~, 

(4.1) \ij_i-i k \\/\ij-i k _i\<N. 
For i, j,k,l > 1, say that 

(4.2) (i,j)~(k,l) 

if 

\(i Aj) - (k Al)\V \(iV j) - (kV l)\ < N . 

Say that a (2m)-tuple (i\, ■ ■ ■ ,12m) is "-/V-pair matched" if there exists a 
pairing tt of {1, ... , 2m} such that 

(ij-i,ij) ~ (V(j-i), V(,-_i)+i) , for j = 1, . . . , 2?n , 

where 7r(0) := 7r(2m). We shall suppress the "iV" in A^-Catalan and iV-pair 
matched if the A^ of interest is clear from the context. 
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Lemma 4.1. Fix N,m,n > 1 and a £ NC2(2m). Let V\, . . . , V m +i denote 
the blocks of the Kreweras complement of a. Write 

K = K, •••,<}, « = l,...,m + l, 

where 

(4.3) <<•••<<. 

TTien i := (ii, . . . , ?2m) £ {1) • • • ,n} 2m is a N -Catalan tuple corresponding 
to a if and only if there exists k := (k±, . . . , fem) £ S(a, N) where 

{Is 
k £ {-N, ..., N} 2m : ^ K? = 0, s = 1, . . . , Jfe + 1 

and a 0- Catalan tuple j := (ji, . . . , j2m) G {1> • • • > n} 2m suc/i that 

X 

(4.5) £„« = j v u + y~] k v n , x = 1, . . . , l u , u = 1, . . . , m + 1 . 

w=l 

Furthermore, the j and k satisfying (4.5) are unique. 

Proof. It is clear that i := (i\, . . . ,i2m) £ N 2m is a 0-Catalan tuple corre- 
sponding to a if and only if 

i v u = . . . = i v y u = 1, . . . , m + 1 . 

In view of the ordering (4.3), similar reasoning as that leading to the above 
equivalence will yield that i is a iV-Catalan tuple corresponding to a if and 
only if 

(4.6) U„« — i„" V i„" — i„" V ... V i,,™ — i,,« < iV, u = 1, . . . , m + 1 ■ 
Now, suppose that i is a iV-Catalan tuple corresponding to cr. Define 

fc„« := i„u +i - i ua , u; = 1, . . . , l u , u = 1, . . . , m + 1 , 

where vf +1 := vf for u = 1, . . . , m + 1. It is easy to see because of (4.6) 
that fe := (fci, . . . , A:2 m ) thus defined, belongs to <S(cr, iV). Define 

i«« := V , w = 1, . . . , l u , u = 1, . . . , m + 1 . 

w lu 

Then, clearly (4.5) holds, and from the equivalence mentioned at the be- 
ginning of this proof, it is easy to see that j := (ji, . . . , j2m) is a 0-Catalan 
word corresponding to a. This completes the proof of the "only if" part. 

For the "if" part, let j and k be given, and define i by (4.5). Then, (4.6) 
is immediate, and by the equivalence mentioned just above that equation, 
it follows that i is a iV-Catalan tuple. 

Finally, for the uniqueness, assume that i is given. If j and k satisfy (4.5), 
then from the fact that 

w=l 
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it follows that 



' u 



This, along with the fact that j is O-Catalan, specifies j. Then, (4.5) deter- 
mines k. This completes the proof. □ 

Lemma 4.2. Let a £ NC2(2m), j be a O-Catalan tuple corresponding to a , 
k £ S(o~,N) and i be given by (4.5). Assume furthermore that 

(4.7) ^{distinct numbers in (ji, . . . , j2m)} = m + 1, 
and that 

(4.8) min{|j u - j v \ : 1 < u, v < 2m, j u ^ j v } > AmN . 
Then, for all (u, v) £ a , 

E {Xi u _ lt i u Xi v _ lt i v ) = R(k u , k v ) . 
Proof. We start by showing that 

(4.9) i u - i v -i = k v 

(4.10) and i v — i u -\ = k v 



^U > 



\, v . 



Assume without loss of generality that u < v. Therefore, u and v — 1 belong 
to the same block in K(a), and furthermore, the block containing them is 
a subset of {u, u + 1, . . . ,v — 1}. Thus, (4.9) follows from (4.5). We show 

(4.10) separately for the cases u > 2 and u = 1. If u > 2, then u — 1 and v 
are in the same block of K(o~), and furthermore that block does not intersect 
with {u, u + 1, . . . ,v — 1}. This shows (4.10), once again with the help of 
(4.5). If u = 1, then 2m and v are in the same block of K{a). Obviously, 
2m has to be the last member of its block. Since (l,f) £ o~, it follows that 
v is the first member of the block containing itself and 2m, showing that 

I'V = *2m i k v = Iq -\- k v = % u —l t k v . 

This complete the proof of (4.10) 
Our next aim is to show that either 

(4.11) % u V i v -\ < i u -i A i v , 
or 

(4.12) i u A i v -\ > i u ~\ V i v , 



holds. Since v — 1 and v belong to distinct blocks of K(a), (4.7) implies 
that 

jv / jv-l ■ 
This, in conjunction with (4.8) establishes that 

\jv ~ jv-i\ > 4mN . 

In view of (4.5), it follows that 

\hi ~ 3v\ V \i u -i -j v \ V \i v -i - j v -i\ V \i u - j v -i\ < 2mN . 
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If j„ > j v -i, then in view of the above two inequalities, it is easy to see that 

i u -i A iv > Jv ~ 2miV > ,/„_! + 2mN > z„_i V i u , 

showing that (4.11) holds. Similarly, if j v < j v -i, then (4.12) holds. 

Finally to see the claim of the lemma, assume that (4.12) holds. Then, 
by stationarity, (4.9) and (4.10), it follows that 

E [Xi u _ lt i u Xi v _ lt i v ) = R{—k v ,—k u ) 
= K{k u , k v ) , 

the second equality following from (2.4). It is easy to see that when (4.11) 
holds, then the claim also holds. This completes the proof. □ 

Lemma 4.3. Fix N > 0, m > 1. In what follows, "pair matched" and 
"Catalan" mean "N-pair matched" and "N -Catalan" respectively. 

(a) I/tt G P(2m) \ NC 2 (2m), then 

lim n~^- m+ >#{pair matched tuples in {1, . . . ,n} m corresponding to tt} 

= 0. 

(b) IfaeNC 2 (2m), then 

lim n~' m+ ># {pair matched tuples in {1, ... ,n} m corresponding to a} 

n— >oo L 

(4.13) \ {Catalan tuples in {1, . . . ,n} m corresponding to a} = . 

(c) Ifae NC 2 (2m) and vr G P(2m) \ {a}, then 

lim n~^ m+ '# {Catalan tuples in {1, . . . , n} m corresponding to a} 



(~]{pair matched tuples in {1, ... ,n} m corresponding to rr} 







Proof. We shall prove the claims only for N = 0. The generalization to the 
case when N ^ is trivial, and follows, for example, by arguments similar 
to the ones that allow replacement of u\ and vi by zero in (3.18) and (3.19) 
respectively. 

Proof of (a). Fix ir G P{2m) \ NC2(2m). Let A = a\ . . . a2 m be a word of 
length 2m such that for j < k, aj = at if and only if (J, k) G n. That is, A 
is a pair matched word. It is easy to see that any 0-pair matched tuple is 
actually an offspring word (considering the entries to be letters) of A. Since 
A is not a Catalan word, by Lemma 3.1 (a), it follows that for all offspring 
word B of A, 

(4.14) #B<m. 

This completes the proof of (a). □ 
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Proof of (b). Let A = a\ . . . a2m be a word of length 2m such that for j < k, 
a j = Ofc if an d only if (j, k) £ a. It is easy to see that a 0-pair matched tuple 
which is not a O-Catalan tuple generates an offspring word B = b\ . . . 62m 
such that at least one of (3.3) or (3.4) is violated for some (j, k) € a. 
Therefore, by Lemma 3.1 (c), (4.14) follows, thus proving (b). □ 

Proof of (c). Once again, we prove this for N = 0. If it ^ NC2(2m), then 
the claim follows by part (a) which has already been proved. So assume that 
7r G NC2(2m). Let A\ and A2 be Catalan words of length 2m corresponding 
to a and ir respectively, as in the proof of (a). Since ir 7^ a, A\ and A2 are 
distinct, that is neither of them can be obtained from the other by relabeling 
letters. Lemma 3.1 (b) implies that if B is an offspring of both A\ and A2, 
then (4.14) holds, thereby establishing (c). □ 

Since all the claims have been established, this completes the proof of the 
lemma. □ 

Remark 2. The number on the left hand side of (4.13) is not necessar- 
ily zero. For example, (1,3,2,1,3,2) is 1-pair matched but not 1-Catalan 
corresponding to {(1, 6), (2, 5), (3,4)}. 

Proof of Theorem 2. 1 . Our first task is to show the existence of a probability 
measure \i whose odd moments are zero and the 2m-th moment is /32m for 
all m > 1. To that end, define the expected ESD /}„, of A n j ^fri as 



1 " 

(B):=-^F(V^e5),n>l, 



n 



for all Borel sets B. Clearly, 

I ' x m (i n {dx) = n-W^E [Tt{A™)\ , m,n > 1 , 

which is zero if m is odd. If it can be shown that for m > 1, 
(4.15) lim n~( m+1 )£ [Tr(A 2 n m )] = (3 2m , 



l (i n {dx) = n-W^E [Tt{A™)\ ,m,n> 



then existence of /i will follow. In addition, the above is also a significant 
step in proving that /x n converges in probability to fj,. We shall come to that 
a moment later. Before that let us quickly dispose off the issue of uniqueness 
of fj>. It will be shown in Section 6 that /?2m. are the moments of a compactly 
supported probability measure which automatically ensures uniqueness; see 
Remark 3. However, for the sake of completeness, we provide a quick proof 
of uniqueness via Carleman's condition. In view of Carleman's condition, it 
suffices to show that 



(4.16) J2 A 



l/2m 
^2r 
m=l 



2m -°° 
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By (2.7), it follows that 

oo oo 

m=l m=l 

1 oo 

> Lr-i/2 J- (m!)- 1/2m 

m=l 
_. oo 

m=l 

the inequality in the last line following from the fact that ml < m 2m for 
all m > 1. This establishes (4.16). Consequently, there is at most one 
measure \x whose odd moments vanish and the 2m-th moment is /32m- Thus, 
to complete the proof, we need to show (4.15) and that 

(4.17) lim n- 2 ( m+1 )Var [Tr(yl 2m )l = . 

n— >oo L J 

We now proceed towards showing (4.15). Recall that 



and therefore 

E[Tr(A*r)]= J2 E " 

ie{l,...,n} 2m 

where 

2m 

E { :=E Y[X i:j _ uij 

3=1 

the convention being that for i = (i\, . . . ,i2m) £ r L 2m , iq := %2m- Recall the 
definition of S(a,N) from (4.4) for a G NC 2 (2m). Set 

& ): = E E II R{K,K),m,N>\. 

<r£NC 2 (2m) keS(a,N) (u,v)&a 

In view of (2.7), it follows that 



i e 



72m 



lim /?£? = /3 2m . 

N— >oo 



Fix e > 0. Let TV" be such that 



\R(u,v)\ < e 



whenever u f\v > N, and 

(4.18) 



P2m ~ Pirn 



<£. 



Recall the definitions of "~", iV-pair matched and iV-Catalan words from 
(4.2) and the text following it. Let PM denote set of the iV-pair matched 



20 A. CHAKRABARTY, R. S. HAZRA, AND D. SARKAR 

tuples in {1, . . . , n} 2m , and let NPM denote the corresponding set for the 
ones that are not TV-pair matched. Write 

E[Tr(A 2 n m )} = £ Bk + J2 E < 

iePM i£NPM 

=: T l +T 2 . 

We shall apply Wick's formula for estimating E{. Fix i E NPM. Recalling 
that P{2m) is the set of all pair partitions of {1, ... , 2m}, given ir E P(2m) 
there exists (u, v) E tt such that 

(iu-i,i u ) 7° (iv-i,iv) ■ 
Therefore, every tt E P{2m) can be written as 

(4.19) vr = {«,<),...,«,<)}, 
where 

(4.20) (iuj-uiiq)'^ (ivj-i,^)- 
Notice that 

(4.21) E[X itj X kt i]=R((iJ)*(k,l)), 

where * is as defined in (3.20). By Wick's formula, it follows that 

m 

' — * ■*■■*■ I L j j j j J 

7r6P(2m)j=2 

m 

7r6P(2m) j=2 

the equality following by (4.21). Therefore, 

TO— 1 

it 2 i < ^ £ #KK«)ni%^i)i. 

(4.22) < eC{m)n m+l R rn - 1 , 

where V^, and C(m) are as in (3.23) and -R is as in (2.3), the last inequality 
following by Lemma 3.7. Thus, 

(4.23) limsupn-( m+1 )|T 2 | < eCim)^" 1 . 

n—*oo 

We now work with T\. For a E iVC2(2m), let 

CT(a) := {Catalan tuples in {1, ... , n} 2m corresponding to a} , 



CT'(a) := CT(a) \ I |J {Pair matched tuples in {1, . . . , n} 5 

\TTeP(2m)\{cr} 



corresponding to tt} , 
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and 

NCT := PM \ i [J CT'(a) J . 

\aeNC 2 (2m) / 

Clearly, 

PM\( (J CT(a)) 

\aeNC 2 (2m) / 

C M {Pair matched tuples corresponding to ir} 

\ir£P(2m)\NC 2 (2m) / 

U I J {Pair matched tuples corresponding to a} 

\a£NC 2 {2m) 

\{Catalan tuples corresponding to a} 

By Lemma 4.3 (a) and (b) respectively, the cardinality of the first and the 
second set on the right hand side is o(n m+1 ). Also, 



( U CT (v))\( U CT\a)) 

\o-6iVC , 2(2m) / \aeNC 2 {2m) / 

|J [CT(a) \ CT'(a)\ . 

a&NC 2 (2m) V / 



7&NC 2 (2m) 

By Lemma 4.3 (c), it follows that the cardinality of the set on the right hand 
side is o{n m+1 ). All the above put together imply that 

(4.24) lim n~( m+1) #NCT = . 

n— >oo 

Clearly, by definition, if a, a 1 G NC2(2m) and a ^ a' , then 

CT'(a) n CT'(a') = <f) . 
Therefore, 

ti = E E ^+ E * 

a&NC 2 (2rn)i&CT'{a) ieNCT 

=■ T n +T 12 . 
By Wick's formula, it follows that 

\Ei\ < (2m)! for all i G N 2m , 
which in view of (4.24), shows that 

(4.25) lim n^ m+1 ^T 12 = . 



22 



A. CHAKRABARTY, R. S. HAZRA, AND D. SARKAR 



Fix a G NC2{2m) and i G CT'(a). By definition of CT'(cr), it is easy to see 
that every n G P{2m) \ {a} can be written as (4.19) such that (4.20) holds. 
Therefore, 



E t - J] E[X iu _ 1>iu X iv . 1>iv ] 

m 

< £ Yl n\R((hi3-hhq)*(ivj-hh,j)) 

■K<=P(2m)\{a} j=2 



By Lemma 4.1, there exist a 0-Catalan tuple j and a k G <S(<7, iV) satisfying 
(4.5). Fix fc G 5(<7, N), n > 4mN and define the sets 



Bi(k) 



B 2 (k) 



i G l? m : (4.5) holds for some 0-Catalan tuple 
j G {1, . . . , n} m corresponding to a > , 



ie{l,...,n} m : (4.5) holds for some 0-Catalan tuple 

j G {1, . . . , n} 2m corresponding to a > , 
B 3 (k) \ i G I? m : (4.5) holds for some 0-Catalan tuple 

j G {AmN + 1, . . . , n — 4miV} m corresponding to a > , 
Bt±(k) = < i G CT (a) : (4.5) holds for some 0-Catalan tuple 

j G {1, . . . , n} m corresponding to a > , 
B 5 (k) = ii G CT'(o-) : (4.5), (4.7) and (4.8) hold for some 0-Catalan 

tuple j G {1, . . . , n} m corresponding to a > . 

Close inspection reveals that 

B 3 (k) C B 2 (k) C B^k) , 
and 



m+l 



#Bi(fc) = n 

#B 3 (k) = (n - 8mN) 



m+l 



Therefore, 

(4.27) 



lim n-( m+ V#B 2 (k) = 1. 
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Lemma 4.1 asserts that 

CT(a) = |J B 2 (k), 

k€S(a,N) 

(4.28) CT'(a) = \J B 4 (k) . 

keS((T,N) 

An outcome of the above is that 

(4.29) B 2 (k)\B 4 (k)cCT(a)\CT\a). 
By Lemma 4.3 (c), it follows that 

lim n"( m+1 )# \CT{a) \ CT'(g)} = , 

which along with (4.27) and (4.29) show that 

lim n^ m+1 ^#B 4 (k) = 1. 

Elementary combinatorics shows that 

(4.30) #lB 4 (k)\B 5 (k)]=o(n m+1 ). 
Therefore, 

(4.31) lim n-( m+1 )#5 5 (£;) = 1 . 

n— >oo 

By Lemma 4.2, it follows that 

(n,i))Scr (u,v)&a 

Similar arguments as those leading to (4.22), in view of (4.26) and the fact 
that the family of sets (B^(k) : k £ S(a,N)) are disjoint, now establish 

E E E *- E #^(£0 [J fl (^^) <eC(m)n m+1 R m ~\ 

k£S(cr,N)i£B 5 (k) k&S{a,N) (u,v)£cr 

where C(m) and R are as in (3.23) and (2.3) respectively. From here, the 
fact that S(a,N) is a finite set, and that (4.31) holds for all k £ S(a,N), 
imply that 



lim sup 



n 



-(m+l) 



E E E i~ E ^"'^ 



< 



keS(a,N) i&B 5 (k) keS(a,N) (u,v)&<7 

Equation (4.30) allows us to replace B§{k) by B±(k) in the above equation 
which along with (4.28) and the observation that the sets on the right hand 
side are pairwise disjoint, lead us to 



(4.32) 



lim sup 



n 



-(m+l) 



e e *- e n R ( k ^ k - 

ieCT'(cr) k£S(a,N) (u,v)€e 
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m—1 



< eC(m)R 
Adding over a S NC2(2m) yields that 



lim sup 

n— »oo 



II 



-(m+l)rp 



R (N) 

~ P2m 



< eC{m)R m - 1 #NC 2 {2m) . 



Since e is arbitrary, the above in view of (4.18), (4.23) and (4.25) complete 
the proof of (4.15). 

To complete the proof, (4.17) needs to be shown, or equivalently, in view 
of (4.15), 



(4.33) 


lim 

n— ¥<x 




m+l)^ 


~{MAi m )f 


— P2m ■ 


Notice that 






E 


~{Tr(A 2 n m )} 2 \ 


= £ * 

ie{l,...,n} 4m 


Am 
J =1 


where 






t(. 


.-{ 


J-1 

2m, 

4m, 


j = 
j = 


{l,...,4m}\{l,2m + l} 

1, 

2m+ 1. 



Denote for i € {1, . . . , n} 4m , 

E; 



E 



4m 

11^70')^ 



Once again, the above expectation can be computed via Wick's formula. 
Therefore, a similar combinatorial analysis as that for the expected trace 
goes through, except that now offspring words are replaced by compound 
offspring words. A sketch of the proof is given below. 

Fix N > 1, and say that % G {1, . . . , n} 4m is iV-pair matched if there exists 
a pairing tr of {1, ... , 4m} such that 



\ij A i l{j) - * ff(j -) A i 7(7r(i)) | V \ij V i l{j) - i nU) V * 7 ( ff (j)) | < N, j 



Am . 



Let PM and NPM denote the sets of iV-pair matched tuples and non-iV- 
pair matched tuples respectively, in {1, . . . , n} 4m with this new definition of 
"pair matched". By arguments similar to those leading to (4.23), Lemma 
3.4 now playing the role of Lemma 3.3, it follows that 



(4.34) 



lim lim sup n 



-(m+1) 



ieNPM 



0. 
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2 r > 



Say that i G {l,...,ra} 4m is Catalan (o"i,a" 2 ) for some 0"i,cr 2 G iVC 2 (2m) 
if (ii, - - - , 12m) is iV-Catalan with respect to a±, and (i 2m +i, ■ ■ ■ , Hm) is N- 
Catalan with respect to <r 2 . By Lemma 3.5, it follows that 



lim n~( m+1 )# I PMA | J CT(cti,ct 2 ) 

,cri,<72e-/VC2(2m) 



0. 



where 



CT(cri,o- 2 ) := {i G {l,...,n} 4m : i is Catalan^, cr 2 )} , ffi,<72 G iVC 2 (2m) . 

By standard combinatorial arguments, it follows from the above equation 
that 



(4.35) 



lim n~( m+1 ) 



x> 



- E E ^ 

iePM a 1 ,cr 2 eNC2(2m) ieCT(<ri,<r 2 ) 



0. 



The arguments leading to (4.32), once again with the aid of Lemma 4.2, 
imply that for o"i, cr 2 G iVC 2 (2m), 



(4.36) 



lim lim sup 

N—too n _>oo 



■». 



-(m+1) 



£ 



£; 



i6CT(CTi,<r 2 ) 



e n *&.*& 



n fl(^,^ 



. (u,«)eo-2 



Clearly, 




E EII *(«)] I n *(&*& 

0-1,0-2 6 7VC 2 (2m) k*eS(<Ti,N) \(u,v)eai J \(u,«)6ff2 

= n( e En *(*u 

»=1 \o-ieAfC2(2m) WeSiiTuN) (u,v)<Scri 

= hm ■ 

This, in view of (4.34) to (4.36), establishes (4.33), and thus completes the 
proof. □ 

5. The linear process 

In this section, we study the ESD of a random matrix whose entries are 
generated from a linear process with independent random variables as the 
input sequence. In particular, let {eij : i,j G Z} be independent, mean zero, 
variance one random variables which satisfy the Pastur condition 



(5.1) 



1 n 
lim -=■ V Mef ,l(k j | > ey/n)] =0 for all e > . 



»,i=i 
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Let {ck,i : k, I G Z} be a collection of deterministic real numbers such that 
(5.2) 0< Y^ \°k,i\ <oo, 



kl&L 



and 

(5.3) 

Define 

(5.4) 



Ck,i = Q,k, k,l ez. 

^i,j := Z^ Ck,l^i-k,j-l, hj G 2, 



where the sum on the right hand side converges in 1? because Ck,i are square 
summable, which is a consequence of (5.2). While the family of random 
variables {Zij : i,j € Z} need not be stationary because the distributions 
of €i j are not necessarily identical, it is easy to see that 

E{Zij) = 0,i,j€Z, 

(5.5) E(ZijZi_ u j +v ) = } c k jc k - u j +v =: R(u,v), i,j,u,v e Z. 

fcJez 

Define the n x n symmetric random matrix ^4 n , and /U n , the ESD of A n /^/n 
by (1.1) and (1.2) respectively. The assumption (5.2) ensures that 

2 



£ix^)i<£x: 



fcez Ze2 



l Cfc .'l z^z_> ICfc-M.Z+^l 



u6/u6i 



zZ l Cfc " 



fc,«e2 



< oo 



Therefore, we can and do define (52m by (2.8). Let \x be the unique proba- 
bility measure whose odd moments are all zero, and for m > 1, the 2m-th 
moment equals ftm. 

The content of this section is the following result. 

Theorem 5.1. Under assumptions (5.1) to (5.3), [i n converges weakly in 
probability to //. 

Proof. We split the proof into two parts. For a finite linear process, we show 
that the Stieltjes transform of the ESD of a matrix made up of Gaussian 
random variables and another with general entries satisfying (5.1) are close 
to each other using Lindeberg type argument developed in Chatterjee (2005). 
For the second part, we show that the Levy distance between the truncated 
linear process and the original process goes to zero as the truncation level 
goes to infinity. 
Fix m > 1 and let 

m 
r (m) 



Z iJ = /Z c k,l£i-k,j-l for i,j > 1 . 



kl=—m 



Define 



4 m) :=((47))nxn,n>l 
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We next define a similar random matrix model, but with Gaussian entries. 
Let (Gij : i,j G Z) be i.i.d. standard Gaussian, and set 

m 

Y ij = E c k,iGi-k,j-i for i,j > 1. 

k,l=— m 

Denote 



4 m) := (0#°))»x», » ^ ! 



"*ij 



Assumption (5.3) ensures that the matrices An and B n are symmetric. 
By the Lindeberg type argument of "replacing eij by Gij one at a time" , it 
can be shown using the Pastur condition (5.1) that 



(5.6) - 
n 



A {m)\ 1 \ // R (m)\ 



Tr *I„ - -^ - Tr z/ n 



n I \ \ wn 



Ao. 



as n — 7- oo, for all z in the complex plane with non-zero imaginary part. The 
arguments for above are very similar to those in Subsections 2.3 and 2.4 of 
Chatterjee (2005) and hence are omitted. 
It is easy to see that 



R^(n,v):=E(Y^Ytl +v )= £ c k , lCk _ u , l+v , u, v G 



k,l=— m 

and as a trivial consequence of (5.3), it follows that 

R(™)(u,v) = R {m \v,u) . 

Clearly, R^ m >{u,v) is non-zero for only finitely many u,v, and hence Theo- 
rem 2.1 applies, with a scaling because R( m '(0, 0) need not be one. By that 
result, it follows that for fixed m, as n — > oo, the ESD of B n / ' \fn converges 
weakly in probability to the probability measure \j> m ' whose odd moments 
are all zero and for I > 1, the 2l-th moment is p^ 1 defined by 

&■■= e e n R {m \k u ,k v ),i>i, 

a£NC 2 (2l) k£S(a) {u,v)ea 

with S(a) being as in (2.6). This, along with (5.6), implies that as n — > oo, 

_1 \ „ r 1 

^ m \dx), zeC\R. 

z — X 




Recalling from (2.9) the definition of L, a restatement of the above is that 

(5.7) i(/ii m U W )A0, 

as n — > oo, where /i„ denotes the ESD of An / y/n. Notice that 
lim R {m) (u, v) = R(u, v) , u, v G Z . 
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By using (5.2) to interchange limit and sum, it follows that 

lim P { ™ ] =(3 2U 1>1. 

m—>oo M 

Therefore, 

(5.8) lim L(V m \//) =0. 

In view of (5.7) and (5.8), to complete the proof of the result, it suffices to 
show that 



(5.9) lim lim sup E L 3 [n { ™>,n 



m— >oo 



0. 



recalling that fj, n is the ESD of A n j ' \fn. 

To that end, we shall use the fact that for n x n (deterministic) symmetric 
matrices C and D with ESD vq and U£> respectively, 

L\v c ,v D )<^Tr{(C-Df) , 

which is a consequence of the Hoffman- Wielandt inequality; see Corollary 
A. 41, page 502 in Bai and Silverstein (2010). Using this inequality, it is 
immediate that 

= 2^ c k,i ■ 

k,leZ:\k\y\l\>m 

The assumption (5.2) of course ensures that {ck,i\ is square summable, and 
thus establishes (5.9). Combining this with (5.7) and (5.8) completes the 
proof. □ 

6. STIELTJES TRANSFORM 

In this section a characterization of the Stieltjes transform of fi, the LSD 
in Theorems 2.1 and 5.1, is given via a functional equation. As the reader 
may have already noticed, in both the above results, [i is defined via the 
correlations R(u,v) which is as in (2.1) or (5.5). For this section, let R(-,-) 
be the correlations of a weakly stationary mean zero variance one process 
(Yij : i,j G Z), that is, 

E(Y id ) = 0,i,jeZ, 
E(Y? d ) = l,i,j€Z, 

E(YijYi- Utj+v ) =: R(u,v), i,j,u,v GZ. 

As in Section 4, we assume (2.2) and (2.3). As before, let [i be the unique 
even probability measure whose 2m-th moment equals /?2m which is as de- 
fined in (2.8). Recall that the Stieltjes transform of the probability measure 
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fj, on R is denoted by, 

G( z ) = / n(dx),z G C 

JrZ-x 

The main result of this section is Theorem 6.1 below. 

Let the Fourier transform of covariance function {R(k,l)}k i<^z be given 

by 

f(x,y) = J2 R(k,l)exp(2m(kx + ly)) for (x,y) G [0,1] x [0,1]. 
k,iez 

Note that by (2.2), it follows that f(x,y) is a real, symmetric function. For 
stating the main result, we need the following proposition. 

Proposition 6.1. Suppose that T~L\ and T-L 2 are functions from C x [0, 1] to 
C satisfying the following for i = 1,2: 

(1) for all x G [0, 1] and z £ C, 

(6.1) zHi{z,x) = l+Hi(z,x) Ui(z,y)f(x,y)dy, 

Jo 

(2) there exists a neighborhood N{ (independent of x) of infinity such 
that for all x G [0, 1], %i(-,x) is analytic on N{, 

(3) for allxe [0,1], 

(6.2) lim zHi(z,x) = 1, 

z— >oo 

(4) and 

(6.3) U(-z, x) = -H(z, x),z€C,x€[0,l]. 

Then 

Hi = H 2 on N 1 r\N 2 . 

Proof. By the assumption of analyticity on N\ and (6.3), for all x G [0, 1] 
and k > 1, there exist H2k(x) G C such that 

oo 

Ui{z,x) =J2 H 2k(x)z~( 2k+1 \ zeN u xe [0, l] . 

fc=0 

The condition (6.2) implies that 

(6.4) H Q (x) = l. 

By comparing the power series expansion of the LHS and the RHS of (6.1), 
one will arrive at the recursion 

m „i 

(6.5) H 2m (x) = y2H 2{m _ k) (x) f(x,y)H 2{k _ l) (y)dy. 

fc=i Jo 

Clearly, a power series expansion of H2 will also satisfy (6.4) and (6.5), and 
therefore they have to match with that of T~L\. This completes the proof. □ 

The following is the main result. 
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Theorem 6.1. There exists a function % satisfying the assumptions of the 
Proposition 6. 1 . The Stieltjes transform Q of the LSD [i is given by 

-l 



9{z) 



W(z,x)dx 
o 



,zeC. 



Before proceeding to the proof, we introduce some notations which will be 
used for the same. Fix a G NC2(2m), and denote its Kreweras complement 
by (Vi, . . . , V m +i). Although the Kreweras complement is a partition of 
{1, . . . , 2m}, for the ease of notation, V\, . . . , V m +\ will be thought of as 
subsets of {1, ... ,2m}, that is, the overline will be suppressed. In order 
to ensure uniqueness in the notation, we impose the requirement that the 
blocks V\, . . . , V m +\ are ordered in the following way. Ifl<i<j<m + 1, 
then the maximal element of Vi is strictly less than that of Vj. Let T a be 
the unique function from {1, . . . , 2m} to {1, . . . , m + 1} satisfying 

* G V T*(i)> I <i <2m. 
For example, if 

a := {(1,4), (2, 3), (5, 6)}, 

then %(1) = 2,%(2) = 1,%(3) = 2,T CT (4) = 4,T CT (5) = 3,7^(6) = 4. 
Define the function L a from M m+1 to M by 

L a {x):= H f(x Ta{u) ,x Ta{v) ),x€W m+l . 

(u,v)£a 



Observe that 






(6.6) 




%(2m 


Finally, set 






K(y) ■ 


Jo 


.. / L ff (xi 

Jo 



,x m ,y)dx m ...dxi, yeR. 

The following lemma shows how the moments of the LSD \x are related 
to this function. 

Lemma 6.1. Let (52m be as in (2.8). Then 

ft™ = ^ / h a (y)dy, m > 1 . 

a&NC 2 (2m) ° 

Proof. Fix m > 1 and a G NC%{^vn). All that needs to be shown is 



(6.7) V R(K,k v ) = f L c 

, 7T, w t ifo,il m+1 



keS(er) (u,v)e<r 



(x) dx , 



where S(a) is as in (2.6). Observe that 

2m m+1 

Yl [ k ^ x T a (u) + Kx Ta{v) ] = Y k u x %(u) =^2x t ^2kj, 
(u,v)€ff «=1 1=1 jeVi 
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and hence 



m + l 



m + l 



L a (x) dx 



m+l 



I 

y, ( n R ( k ^kv)) [ 



^ exp 27ri ^ x z ^ kj R(k u , k v ) 



da; 



m+l 



exp 2-7ri N^ x/ N^ fcj dx 



= Yi n ^a), 

fceS(o-) (u,i))e<r 

the interchange of sum and integral in the second last line being justified by 
assumption (2.3). This completes the proof. □ 

The next two lemmas give a recursive relation for the function h a . 

Lemma 6.2. Assume that a G NC2(2m) can be written as 

(6.8) a = a x U a 2 , 

where u\ £ NCi^lK) for some 1 < k < m — 1, and 02 is a non-crossing pair 
partition of {2k + 1, . . . , 2m}. Viewing o~2 as an element of NC2(2m — 2k) 
by the obvious relabeling of 2k + 1, ... , 2m to 1, . . . , 2m — 2k respectively, it 
is true that 

My) = M(j/)M(y)> v G M - 

Proof. It is easy to see from (6.6) and (6.8) that 

(6.9) %(j) e{l,...,k,m + 1}, for 1 < j < 2k , 
and 

(6.10) %(j) e{k + l,...,m + l}, for 2k + 1 < j < 2m . 

Write 

L a (x) = 

n /(^oo^m) j i n / i x %(u),x%(v)) 

(u,v)£a:u,v<2k J \(u,v)£a:u,v>2k 

By (6.9) and (6.10), it follows that 
h ff (Xm+l) 

= [j ■j f{x%(u),XT*(v))dx k ...dxi\ 

\ (u,v)da:u,v<2k J 

(/■■■/ ^ ( x r CT («)^r CT (i>)) da?m • • • dx k+1 

\ (u,u)S<7:ivu>2fc / 
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This completes the proof. □ 

Lemma 6.3. If 

a = {(1,2m)} Ucti , 
for some non-crossing pair partition o\ of {2, ... , 2m — 1}, then 

h„(z)= K 1 (y)f(y,z)dy, z el, 
Jo 

where, once again, o\ is viewed as an element of NC2(2m — 2). 

Proof. Throughout the proof, o~\ is to be thought of as an element of 
NC 2 {2m-2). Clearly, 

%{l) = m, 
T a {2m) = m + 1 , 

%(j) = % 1 (j-l),2<j<2m-l. 

The above equations imply that 

L a (x) = f(x m ,x m+1 )L ai (xi, . . . ,x m ), x £ R m+1 . 
Thus, 

K(z) = / •••/ f(x m ,z)L ai (xi,...,x m )dx m ...dxi 
Jo Jo 

= / f(x m ,z)l ... L ai (xi,...,x m )dx m -i...dxi\dx m 
l 



which completes the proof. □ 

With the aid of the above lemmas, we proceed to the proof of Theorem 6.1. 
Proof of Theorem 6. 1 . We start with defining the functions 
H Q {x) = 1, H 2 m(x) = ^2 h„(x). 

aGNC 2 {2m) 

Lemma 6.1 implies that 

(6.11) hm= f H 2m {x)dx. 

Jo 

By (2.3), it follows that \f(x,y)\ < R uniformly for all x and y in [0,1]. 
Lemmas 6.2 and 6.3 applied inductively imply that 

sup \K(x)\ <R m ,a£ NC 2 (2m), m > 1 . 

0<x<l 
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Since #NC 2 {2m) < 4 m , it holds that 

sup \H 2m (x)\ <R m 2 2m ,m> 1. 

0<x<l 

Consequently, 

(6.12) sup limsup\H^ m (x)\<2R 1/ ' 2 <oo. 

0<x<l m— >oo 

Therefore, the power series 

n(z,x) = ^^^r 

m=0 

converges on {z G C : \z\ > 2.R 1 ' 2 } for every fixed x G [0,1]. Note that 
this neighborhood around infinity is independent of x £ [0,1]. It is easy 
to see that zH{z,x) has a power series expansion with the leading term as 
1 and hence zH(z,x) -f 1 as \z\ — > oo. It follows from the definition of 
H(z,x) that H(—z,x) = —H(z,x). Recall that the Stieltjes transform Q of 
/j, satisfies 



Q{z) = ^ ftmZ" 



(2m+l) 



m=0 

(with the obvious convention that /3q := 1) which yields, 

g{z) = [ U{z,x)dx. 
Jo 

Equation (6.1) with Tii replaced by % is all that remains to be checked. 

To that end, we derive a recursion for H(z,x) using the properties of 
h a {x). Recall that there is a natural one-one correspondence between 
NC2(2m) and the set of Catalan words of length 2m with the understanding 
that two words will be considered identical if one can be obtained from the 
other by a relabeling of letters. Keeping this correspondence in mind, by 
an abuse of notation, we shall now consider h w (x) for Catalan words w, 
and denote by ./VC2 (2m) the set of Catalan words of length 2m. Note that 
any Catalan word w of length 2m can be written as w = awiaw2, for some 
W! G NC 2 (2k - 2) and w 2 G NC 2 {2m - 2k). So if 



then 



H 2m ,k{x) : = ^2 ^2 h 

w 1 GNC 2 (2k-2) w 2 eNC2(2m-2k) 
in. 

H 2m {x) = ^2H 2m ^ k (x). 

k=l 



awiau)2 \X) 



Notice that 



H 2m ,k( x ) = ^2 K Wl a(x) ^2 h W2 (x) 

W!£NC2(2k-2) w 2 eNC2(2m-2k) 
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= ^2 / if( x ^y) h wi(y)]dy Yl h w 2 {x) 

w 1 eNC 2 (2k-2) J ° w 2 £NC 2 (2m-2k) 

= / [f{x,y)H 2 ( k -i)(y)H2(m-k)(x)]dy, 
Jo 

the equalities in the first two lines following from Lemmas 6.2 and 6.3 re- 
spectively. As a consequence, 

m „! 

H2m(x) = y2H 2 ( m - k )(x) / f(x,y)H 2 (k-i)(y)dy. 
fc=i ^ 

Now by an easy computation it follows that, 

oo 1 1 r 1 

U(z,x) = V H 2m z^ 2m+1 ^ = - + -H(z,x) / f(x,y)H(z,y)dy. 

m=0 Z Z J 

Hence, 

zn(z,x) = l + n(z,x) / f(x,y)n(z,y)dy. 



This completes the proof. □ 

Remark 3. Equations (6.11) and (6.12) imply that 

lim j3l^ m < oo , 
which implies that the probability measure /i is compactly supported. 

7. Special cases and examples 

In this section, we attempt to give a better description of the probability 
measure fi, which appears as the LSD in Theorems 2.1 and 5.1, in some 
special cases. As in Section 6, R(-,-) are the correlations of a weakly sta- 
tionary mean zero variance one process (Yij : i,j E Z). As always, (2.2) 
and (2.3) are assumed, and \i is the unique even probability measure whose 
2m-th moment equals f3 2m which is as defined in (2.8). The first main result 
of this section is the following. 

Theorem 7.1. Assume that 

(7.1) R(u, v) = R(u, 0)R(0, v),u,v€Z. 

Then, the function r(-) defined on [— ir, ir] by 

oo 

r{x):= Y, R{k,0)e~ ikx ,-ir<x<ir, 

k=—oo 

is a well defined function, that is the sum on the right hand side converges 
absolutely, and its range is contained in [0, oo). Furthermore, 

(J, = fJ, r M H a , 
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where \i r denotes the law ofr(U), U being a Uniform (— ir,ir) random vari- 
able, ji s denotes the WSL whose density is given by 



x/4" 



x 



2 



(7.2) fj, 3 (dx) := ^— l(|x| < 2)dx , 

Z7T 

and " M " denotes the free product convolution. 

Remark 4. In Bercovici and Voiculescu (1993), the free multiplicative con- 
volution of two probability measures with possibly unbounded support, at 
least one of which is supported on a subset of [0, oo), has been defined. 
Hence, in the above result, the claim that r(-) is non-negative is needed. 

Proof of Theorem 7.1. In view of assumptions (2.3) and (7.1), it is easy to 
see that 

oo 

(7.3) E \ R ( k '°)\ =R l/2 <oo, 

fc=— oo 

and hence the infinite sum defining r(x) is absolutely convergent. Observing 
that (R(k, 0) : k G Z) is the autocovariance function of the one-dimensional 
process (Y"-fc,o : k G Z), Corollary 4.3.2, page 120 of Brockwell and Davis 
(1991) implies that 

r(x) G [0, oo), — it < x < ir . 
This, in view of (7.3), establishes that the range of r(-) is a compact subset 
of [0,oo). 

For establishing the other claim, we shall use Theorem 14.4 of Nica and Speicher 
(2006) which applies to compactly supported probability measures. From 
that result, it follows that fi r Kl fi s is an even probability measure, and 

/m+l „ 
x 2m (fi r m fx s ){dx) E / x l 3fi r (dx), m > 1 , 

a£NC 2 (2m) j = l J 

where for any a G NC2(2m), If, ■ ■ ■ ,1^+1 denote the block sizes of the 
Kreweras complement of a. It is easy to see from the definition of r(-) that 

r j 

x^ r (dx)= e n^^' )'^ 1 - 

{ki,...,kj)eU :ki+...+kj=0 i=l 

For a G NC2(2m) let the notation for its Kreweras complement be as in 
(2.5), and recall the definition of S(o~) from (2.6). The above two identities 
put together imply that for m > 1, 

/2m 
x 2m ^ r M^ s )(dx) = E E 11^'°) 

aGNC 2 (2m) keS(a) j=l 

E E II R(ku,0)R(k v ,0) 

aeNC 2 (2m) fceS(o-) (u,v)&a 
= ftm 
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x 2r y(dx) , 

the equality in the second last line following from (7.1). This completes the 
proof. □ 

The next result is the other main result of this section. 

Theorem 7.2. // 

(7.5) R(k, 0) = for all k / , 
then fj, = fi s , where [i s is the WSL as defined in (7.2). 

Proof. Clearly, it suffices to show that for all m > 1 and a G NC2(2m), 

(7.6) J2 II R(ku,k v ) = l, 

k&S(a) (u,v)ecr 

where S(a) is as in (2.6). This is because, if the above is established, then 
it will follow that 

fo m = #NC 2 (2m) = J x 2m fi s (dx) . 

In order to show (7.6), fix m > 1 and a G NC2(2m). What we shall show is 
that if k G S(a) is such that 

(7.7) Yl R(k u ,k v )^0, 

(u,v)£a 

then, 

(7.8) k 1 = ... = k 2m = 0. 

Recalling that R(0, 0) = 1, which is a consequence of the assumption that 
the process (Yij : i,j G Z) mentioned at the beginning of this section, has 
variance one, the above will imply (7.6). 

The claim (7.8) is a tautology when m = 1. As the induction hypothesis, 
we assume that for a fixed m > 1 and all a G NC2(2m), (7.7) implies (7.8). 
To complete the induction step, fix a G NC2(2m + 2), and let (7.7) hold. By 
the property of non-crossing pair partition, there exists j G {1, . . . , 2m + 1} 
such that (j,j + 1) G o. Recalling that K(o~), the Kreweras complement of 
a, is the maximal partition a of {1, ... , 2m} such that a Da is a non-crossing 
partition of {1, 1, . . . , 2m, 2m}, it follows that (j) G K{a). Hence, 

kj = 0. 

Since (7.7) holds, it follows that 

#(%,%+!) ^0, 

which along with (7.5) implies that 

k j+1 = . 
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If a denotes the element of NC2(2m) obtained from a by deleting (j,j + 1) 
and the obvious relabeling, then it is easy to see that 

(fa, ... , k 2m ) ■= (fa, . . . , kj-i,kj+2, ■■■ , k 2m +2) e S(a) , 
and 

J] R(k u ,k v )^0. 

(u,v)€a 
By the induction hypothesis, it follows that 

Aii = ... = hm = , 
which establishes the induction step, and thereby completes the proof. □ 

Now, we shall see the relevance of the two main results proved above in 
the light of Theorem 5.1. 

Corollary of Theorem 7.1. Let {e^j : i,j G Z} be as in Section 5; in 
particular, the Pastur condition (5.1) holds. Let {ck ■ k G Z} be a sequence 
of real numbers such that 

oo 

(7.9) Yl N <00 ' 

k=— oo 

and 

oo 

(7.10) Y, °l = l - 

k=—oo 

Set 

Ck,i ■= c k ci, faleZ. 

Define Zij and R(-,-) by (5.4) and (5.5) respectively. Clearly, (5.2) and 
(5.3) hold, and the process (Zij : i,j G Z) is weakly stationary with mean 
zero and variance one. Also, 



R(u, v) = I y^ c k c k -u Y ClCl+v 



c l c l+v c l' 



V k k' J \ I V / 

= R(u,0)R(0,v) , 

the second equality following from (7.10). Let A n and fx n be as in (1.1) 
and (1.2) respectively, that is, the former is the n x n matrix whose (i,j)-th 
entry is ^Aj.ivj'j an d the latter is the ESD of A n j \fri. Let \x T and \i s be as in 
the statement of Theorem 7.1. Then, as a corollary of the result mentioned 
above and Theorem 5.1, it follows that, \i n converges weakly in probability 
to fj, r M fj, s . 
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Corollary of Theorem 7.2. Once again, let {ejj : i,j G Z} be as in 
Section 5 satisfying (5.1). Assume that {c^i : k,l S Z} C R is such that 
(5.2) and (5.3) hold, and furthermore 

oo 

(7.11) 2, c k,iCk',i = 1(A; = k') for all k, k' € Z. 

/ = — oo 

As before, let Zj ,-, ^4 n and /x n be as in (5.4), (1.1) and (1.2) respectively. It 
is easy to see that the conditions imposed above ensure that (Zjj : i,j £ Z) 
is a mean zero variance one weakly stationary process, and that (7.5) holds. 
Then by Theorem 5.1 and Theorem 7.2, it follows that \x n converges weakly 
in probability to p s which is the WSL defined in (7.2). 

We end this section by revisiting Examples 1 to 4 mentioned in Section 
1. 

Example 1. To start with, one needs to argue the existence of a stationary 
centered Gaussian process {Zij : i,j E Z} satisfying 

E[Z 0> oZu iV ]=p\ u \ + \ v \,u,veZ. 
That, however, is obvious from the observation that 

pH+M = f e i{ux+vy) F{dx)F{dy), u,veZ, 

J(-7T,7r] 2 

where F is the spectral measure of the autocovariance function (p' ' : h £ 
Z); see Herglotz theorem (Theorem 4.3.1 in Brockwell and Davis (1991)). 
By Theorem 7.1 and results about the AR(1) process, it follows that p n 
converges in probability to p r K p s , where p r is the law of -r-^ — ~ p f/ 5 , U 
being an Uniform (— vr, n) random variable. 
Example 2. Notice that the (i, j)-th entry of A n is given by 

(N + 1) J2 cmGi-^j-i =: (N + l)Yi,j , 

k,iez 

where c k := (N + l)~ 1 / 2 l(-Af < k < 0). Then (7.9) and (7.10) hold, and 
therefore, the ESD of ((Yij / '\^n)) n xn converges to p r M p s , where p r is the 
law of 

N 

1 + 2(N + l)" 2 J2( N -k + lf cos(kU) , 
fc=i 

U being distributed as Uniform (— 7r,7r). Hence the LSD of A n /sjn is the 
free product convolution of p a with the law of 

N 

N + 1 + 2(N + ly 1 J2(N -k + lf cos(kU) . 
fc=i 

Example 3. By Theorem 7.2, it follows that under the additional assump- 
tion that Yln=i \E(G G n )\ < oo, the LSD of A n /y/n is p s . 
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Example 4. Setting 

(oo oo \ 1/2 

k=— ool=— oo / 

it is easy to see from the Corollary of Theorem 7.2 that the LSD of 
a^ 1 A n j ' yfn is [i s . Therefore, the LSD of A n / ' yfn is \x a given by 



*/4 — X^ I (7^ 

fl s (dx) := 2 — - — l(|x| < 2a)dx, 

lira 

which is a dilation of the WSL. 

8. AN EXTENSION OF THEOREM 7.1 

In this section, we generalize Theorem 7.1 to the case when the covariances 
R(u,v) are not necessarily summable. Suppose that (Zij : i,j G Z) is a 
stationary mean zero variance one Gaussian process. Define A n , fj, n and 
R(-, •) by (1.1), (1-2) and (2.1) respectively. The first assumption is, as 
before, that 

(8.1) R(u,v) = R{u, 0)R(v, 0), u, t> G Z . 

The second assumption, the one that replaces the summability of R(u,v), 
is that the spectral measure of the one dimensional process (^,o : « G Z) 
is absolutely continuous with respect to the Lebesgue measure. This means 
that there exists a non-negative integrable function r on [— ir, ir] such that 

e inx r(x)dx = R(n,0), n£Z. 

Strictly speaking, (2-7r) _1 r(-) is the density of the spectral measure. Since 
R(n, 0) = R(—n,0), it follows that r(-) is symmetric. As in Theorem 7.1, 
denote by \x r the law of r(U) where U is an Uniform(— ir, it) random variable. 
The main result of this section is the following. 

Theorem 8.1. As n — > oo, fi n converges weakly in probability to /i r Kl /j, s , 
where fj, s is the WSL. 

For the proof of Theorem 8.1, we shall need the following lemma, which 
is an observation of independent interest. 

Lemma 8.1. Define 

(8.3) c fc := (2ny 1 I e ikx \frjx) dx , k€Z, 



the integral being defined because yV(-) G L 2 ([—tt,tt]), and real because r(-) 
is symmetric. Then, 

oo 

(8.4) Y, °l < °° ' 
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and thus the sum Yliki^i c k c iGi-k,j-i converges in L? for all i,j, where 
(Gjj : i,j > 1) is a family of i.i.d. standard Gaussian random variables. 
Furthermore, 

(Z iyj :i,jeZ)= ^2 c k ciGi- k ,j-i :i,je 

Proof. The claim (8.4) follows from an application of the Parseval identity. 
Therefore, 

^ CkCiGi-kj-i :i,jeZ 

is clearly a mean zero stationary Gaussian process. To check that it has the 
same finite dimensional distributions as (Zij : i,j eZ), it suffices to verify 
that the correlations match, that is, 

oo oo 

(8.5) R(u,v) = ^2 /^2 c k°k-uCiCi +v , u, v e Z . 

k=— oo l=— oo 

To that end, we start with the observation that 

(8.6) lim V c k e~ ikx = Jrjx), -vr < x < vr , 

7V->oo *-^ 
k=-N 

in L 2 ([— 7r, 7r]), which follows from the fact that the Fourier series of a square 
integrable function converges in the L? norm to that function. Therefore, the 
square of the left hand side converges in L 1 to r(x) as N — > oo. Consequently, 
for fixed n G Z, 

' N s 2 



f" e inx r(x)dx = lim f e mx I V c k e~ lkx ) dx 

N 

= lim 2?r Y] c fc c n _ fc 

k=-N 
oo 

= 2vr 2_/ C k C n _ k 



k=—oo 
oo 



= 27T 2J c kCk-n , 
k=—oo 

the second last equality following from (8.4), and the last equality being an 
outcome of the fact that y / r(-) is symmetric. Comparing this with (8.2), it 
follows that 



R(n,0) = ^2 c kCk-n, n e 



k=— oo 
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Using (8.1), it follows that for all u, v gZ, 

R(u,v) = ^ c kCk~u\ ^2 C l C l-v) 

\k=—oo / \l=—oo / 

= I ^2 c k c k - u ] I ^2 c i+v°i ) > 

\fc=— oo / \l=— oo / 

thereby establishing (8.5). This completes the proof. □ 

Proof of Theorem 8.1. Let (Gij : i,j > 1) be a family of i.i.d. standard 
Gaussian random variables, and let {c k } be as in (8.3). In view of Lemma 
8.1, without loss of generality, we can and do assume that 

Zij = ^2 c kCiGi~k,j-u i,j £%>■ 
k,iez 

Denote 

m m 

Z if := 5Z ^2 c kClGi-k,j-l, i, J G Z, m > 1 , 



PS m \u,v) := e[z$Z 



k=—m l=—m 

, u, v € Z, m > 1 , 



m) y(m) 
"-u,v 



4 m) := ((#))nxn,n,m>l. 



'..; 



If Hn denotes the ESD of A n /y/n, then exactly same arguments as those 
in the proof of Theorem 5.1 show that 



(8.7) lim lim sup £ L 3 /4 )/M = °> 

m->oo n ->oo I \ /J 

where L denotes the Levy distance. Clearly, 

(8.8) lim R ( - m \0,0) = R(0,0) = l. 

m— >oo 

Fix m large enough so that R^ m '(0,0) > 0. By Theorems 2.1 and 7.1, it 
follows that 

(8.9) L^,^®^ Ao, 
as n — > oo, where 

and /i rm denotes the law of r m (U), U being an Uniform (—ir,ir) random 
variable. Notice that in the sum on the right hand side, only finitely many 
terms are non-zero. In view of (8.7) and (8.9), to complete the proof it 
suffices to show that 

(8.10) lim L{fi a ^fir m ,fi a ^fi r ) = 0. 
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To that end, define 



oo / m \ ^ 




r m (x): = J2 R {m) (n,0)e- inx =i £ c k e~ ikx \ , 


— 7T < X < 7T 


7i=— oo \k=-m J 




By (8.6), it follows that 




lim / \r m (x) — r(x)\ dx = , 

m— 5>oo / 




which along with (8.8) ensures that 




w 

/j, rm — > ^r, as m — > oo . 





Define the map y/~- from the space of non-negative probability measures 
to that of symmetric probability measures as follows. If a non-negative 
probability measures v is the law of a random variable V, then ^/v is the 
law of eVV, where e takes values +1 and —1 each with probability 1/2 
independently of V. Let - 2 denote the inverse of yf~- . By Corollary 6.7 of 
Bercovici and Voiculescu (1993), it follows that 

ry ry qjn ry r* 

Vs E /V m — ► /u s ^ M r , asra^oo. 

Lemma 8 of Arizmendi and Perez- Abreu (2009) tells us that for a symmetric 
probability measure v\ and a non-negative probability measure v^ such that 
*/i({0})V*/ 2 ({0})<1, 

z/1 IE1 z^ 2 = v/^i ^ ^2 • 
This shows (8.10) which, along with (8.7) and (8.9), completes the proof. □ 

Next, let us see two examples where Theorem 8.1 applies. 

Example 5. Let 

sin u sin v 
R(u,v) := , ii,»eZ, 

U V 

where (sin0)/0 is to be interpreted as 1. Clearly, 

(27T)- 1 / e mx irl(\x\ < l)dx = R(n,0), neZ. 

J ~n 

Therefore, by Theorem 8.1 and the fact that the free product convolution 
of Bernoulli (p) and WSL is same as the classical product convolution of 
yfp times Bernoulli (p) and WSL, it follows that the LSD in this example 
is the law of UW, where W follows WSL, and LT takes the values \pn and 
with probabilities 1/ir and 1 — 1/ir respectively, independently of W in the 
classical sense. Note that for this example, even though 

^2\R(u,v)\ =oo, 

u.v 

the LSD is still compactly supported. 
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Example 6. Define 



r{ x ) := ^-\x\~ 1/2 l{x / 0), -vr < x < vr , 
and 
*(_„,:= (<*)-' />*,*) (<*)-/W>*),.,„ 6 Z. 

The fact that 

/n 
r(x)dx = 2-7T , 

ensures that the assumptions (8.1) and (8.2) hold. Therefore, the LSD 

in this case is \i r M fi s , where fi r is the law of ^-\U\~ 1 ' 2 and U follows 
Uniform(— 7r,7r). The following proposition shows that the fourth moment 
of the LSD is infinite, which means, in particular, that the largest eigenvalue 
of A n /y/n goes to infinity in probability. 

Proposition 8.1. For any non-negative probability measure u and integer 
k> 1, 

(8.11) / x 2k /j, s ®v(dx) = oo, 

Jr 

if and only if 

(8.12) / x k u{dx) = oo. 



Proof. We start with the "if" part, that is, assume (8.12). Let X be a 
random variable whose law is v. For n > 1, let v n denote the law oi X An. 
In what follows, all integrals are on the whole of M. By (7.4), considering 
the element {(1, 2), . . . , (2k - 1, 2k)} of NC 2 (2k), it follows that 

x /j, s ^ v n (dx) > I / xv n (dx) J / x u n (dx) . 

Proposition 4.15 of Bercovici and Voiculescu (1993) implies that fi 2 Kl i/ 2 is 
dominated by fj, 2 Kl v 2 , and hence 






> I / xv n (dx) I / x v n (dx 



k 

k 



I xv(dx) J / 



x v(dx) = oo , 



the limit in the last line following from the monotone convergence theorem. 
This establishes the "if" part. 
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For the "only if" part, assume that 



/ 



x v(dx) < oo . 



Define v n as before for n > 1. For a € NCi(2k\ let /J, . . . ,^ +1 be as in 
(7.4). By the Skorohod embedding and Fatou's lemma, it follows that 

/ x 2k fj, s M u(dx) < liminf / x 2k ^ s Kl v n (dx) 
J ™->°° J 

k+1 ,. 
= liminf y \ x i v n (dx) 

n— >oo * — ' -'--'-/ 

aG7VC 2 (2fc)i=l ,/ 

fe+1 r 

<TeAfC 2 (2fc)i=i^ 
< oo , 

the inequality in the last line following from the observation that l? < k for 
all a £ NC2{2k) and 1 < j < k + 1. This completes the proof of the "only 
if" part. □ 
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